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Abstract 

We examine the NMR experimental data for the normal state of 

YBa2Cu307_,5 within the Fermi-liquid approach. We show that the observed 

temperature dependence of the Knight shift and 17 ( 89 Y) relaxation rate in 

"underdoped" ( T c = 60 K ) compound can be interpreted as resulting from 

a peak in the density of states near the Fermi level. The possible origin of 

such peak is the quasi-one-dimensional van Hove singularity in the fermionic 

spectrum, conjectured in Q. The proposed spectrum allows us to account 

for the qualitatively different NMR data for "underdoped" and "overdoped" ( 

T c = 90 K ) compounds by varying the value of chemical potential only. 
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Last years for the description of properties of high-T c cuprates a series of theoretical 
concepts have been proposed, invoked to explain the high-T c phenomenon itself and the 
anomalies of the normal state of these compounds. We wish to mention here the van-Hove 
singularity (vHs) scenario and nearly antiferromagnetic Fermi liquid (NAFL) H model. 

In this paper we address mainly to the NMR experiments in the normal state of un- 
derdoped compound YBa2Cu306.6 JU1 ( ~ 60 K ), which exhibits the properties fairly 
distinctive from those of overdoped YBa2Cu306.95- The main difference relates to the tem- 
perature dependence of the Knight shift above T c . It has been experimentally observed 
@H, that the Knight shift K in underdoped compound increased with temperature and 
apparently saturated at T ~ 300 K , while it should be constant for ordinary Fermi liquid. 
The 17 and 89 Y relaxation rates (TTi) -1 behaved similarly to the Knight shift while 
the copper relaxation rate ^(TTx) -1 suffered a maximum at T « 140 K. 

It was shown || that the qualitative difference in the relaxation rates at different sites, 
as well as the unusual proportionality of (TTx) -1 to K in the underdoped compound can be 
understood in the frame of semiphenomenological NAFL approach. However, the Knight 
shift was treated as a parameter of the theory, i.e. has not been explained. Another pa- 
rameter, the antiferromagnetic correlation length essentially varied with temperature, which 
clearly contradicts to the inelastic neutron scattering data (INS) ||. Next, this theory can 
not be naturally continuated into the superconducting temperature region. 

The temperature dependence of the Knight shift as well as some of INS data for the 
underdoped materials are usually referred as the opening of spin-gap with lowering of the 
temperature 0. 

In the present paper we explain the observed temperature dependence of the Knight shift 
and O ( Y ) relaxation rate for YBa2Cu30g.6 compound within the Fermi liquid approach. 
Our consideration utilizes the hypothesis of the peak in the density of states (DOS) near 
the Fermi level in this system, which idea is shared by the vHs models |2[. This idea has 
been previously discussed in connection with the thermodynamic measurements in YBaCuO 
compound @. On the other hand, the recent photoemission data witness in favor of flat 
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parts of electron spectrum near the Fermi surface |],|8| in high-T c cuprates, which should 
correspond to the peak in DOS at energies close to the Fermi level. 

Leaving apart the widely discussed question whether it is possible to explain high T c in 
the frame of the peak-in-DOS hypothesis, we confine ourself with the consequences for the 
normal state quantities observable in NMR experiments. 

Our main idea is as follows. The values of the spin part of Knight shift K and the 
reduced nuclear spin-lattice relaxation rate (TTi) -1 (RR) in metals are determined by the 
electrons, whose energies ( measured from the Fermi level ) are of order of temperature. 
If the density of states within this interval is nearly constant one obtains the temperature 
independent values of K and (TTi) -1 . However, if the DOS function possesses an additional 
(small) energy scale near the Fermi level, it might result in the essential variation of these 
quantities with temperature. We found that the experimental data for the Knight shift for 
YBa2Cu306.6 are consistent with a peak in DOS at the energy ~ 40meV. On the other 
hand, the relaxation rate is determined by the integrated square of DOS, in which quantity 
the peak is reproduced as well. So one can conclude that the temperature dependence of 
(TTi)^ 1 should be similar to one of K in accordance with experimental data for 17 and 
89 Y. Discussing the origin of such peak in DOS, we argue that the most likely reason for it 
could be the quasi-one-dimensional vHs conjectured previously by Abrikosov, Campuzano 
and Gofron [l[ on the base of photoemission data. 

Moreover, we show that the proposed DOS function permits us to account for the quali- 
tative difference between the compounds with T c m 60 K and T c « 90 K by the difference in 
the chemical potentials only. 

Knight shift. The spin part of the Knight shift is proportional to the static susceptibility 
of electron system. In the normal state it can be represented through the DOS function p(x) 
in the form ||: 




(1) 
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P( X ) =J2 S ( E k- X )\ U k 



(2) 



k 



here A a (q) is the hyperfine coupling constant, 7 e (7^) is the electron (nucleus) gyromagnetic 



the amplitude of Bloch function on the nucleus ( for s-electrons ) or the coefficient in the 
diagonalizing Hamiltonian ( for multi-band electronic structure ) or the effects of carriers' 
interaction. Our main assumption henceforth is small dependence of ( and therefore p(x) 



The usual assumption is that the density of states is almost constant near the Fermi 
surface. It is valid, e.g., for the Fermi gas, when E\~ = k 2 /(2m) and = 1 so that p(x) = 
p(Ep){l + 0(x/Ep))- Hence, for the typical values of the Fermi energy Ep — p ~ 1 eV the 
temperature corrections to K a are of order of (T/Ep) 2 <C 1. 

If however p(x) has some prominent feature at low energy comparable with temperature, 
then the temperature dependence of Knight shift would not be so trivial. Namely, let us 
assume that p(x) in ([!]) has a peak at x — Ei + p. Then the weight of this peak ~ 
(Tcosh^i^T))- 1 rapidly increases, as the temperature increases from zero to T m O.6E1, 
and starts to decrease slowly upon further temperature raising. Obviously the resulting 
Knight shift would have the similar temperature behavior. 

We present below a procedure which may let one judge about spectral features from 
the Knight shift data. We note that the expression (|l|) can be rewritten in terms of Mellin 
transforms in the following way : 



ratio and p is the chemical potential. The quasiparticle weight \uk\ 2 in (0) includes either 



) on T fnj. 



K(s)=p + {s)<p(l-s) 



(3) 



Here we use the following definitions : 



p + (x) 



tp(x) 



2- 1 cosh" 2 (x/2), 
~T~^ 0°(/ i + x ) + P(p ~ x )) 




(4) 



and 
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f(s)= / dxx s f(x), s = a + i\ (5) 
Jo 

f(x) = (2m)- 1 dsx~ s f(s) (6) 



Of course, trying now to obtain from ([3]) the function p + one meets the difficulties 
common for the badly-defined inverse problems, when (p(l — s) — *■ 0. Therefore we do not 
pretend to extract the exact form of p+(x) and find instead only the regularized solution 
p+(x, 5) of the inverse problem in the form : 

where 5 > and a is chosen from the condition of simultaneous existence of Mellin transforms 
K(a) and tp(l—a). To satisfy this condition, we subtract the constant term ( at T > 300 K ) 
from the observed temperature dependence of Knight shift, using the property J °° dx ip(x) = 
1. We note that the restored function p+(x, 5) is almost insensitive to the choice of a G [0, 1] 
and depends strongly on 5. The values of S = (10~ 2 -r 10~ 4 )|</?(1 — a)\ 2 give a shape of 
p+(x, 5) with a visible details and , on the other hand, free from the giant oscillations. 

Choosing o = 0.35 and 5 = 10~ 3 |<^(1 — o~)| 2 , we find the regularized solution for symmet- 
rical part p + of DOS which is shown in Fig. [I] by the dashed line. The resulting DOS has a 
peak at E\ ~ 40 meV. Obviously the fine features of DOS structure could not be restored 
during the procedure of regularization so this picture delivers only the rough information 
about the position and the general shape of the DOS peak. 

Let us discuss briefly the reasons of such strong variation of DOS function at such a 
small energy scale. First possibility relates to the existence of vHs near the Fermi level. 
This possibility was widely discussed before, mainly for the case of logarithmically divergent 
DOS, pix) ~ ln(A/|x — This kind of vHs can be found in a variety of papers based 

on the tight-binding Hamiltonian and originates from the saddle point of two-dimensional 
fermionic spectrum P,|12||. We wish to stress, however, that the logarithmic singularity is 



too weak to produce the essential temperature variation of the Knight shift (cf . |12| ) . 

At the same time, the recent photoemission data for YBaCuO compounds give a 



strong confidence to the existence of the "extended" saddle point along the T — Y symmetry 
direction. It has been shown [|TJ that the electronic spectrum near the (0, 7r) point of the 2D 
Brillouin zone can be represented in the form 

E k = tk 2 x + Ei, \k y \ < k ps .5r.l.u., 

where E\ was found to be within the experimental resolution of 30 meV. The above kind of 
spectrum results in DOS function of the form 

p(x) = p n + p,(x), (8) 

where p s (x) ~ p n \jtk\j [x — E\) 9(x — Ei), and p n ~ t -1 is the smooth part of DOS. The 
divergency in p(x) is obviously smeared out either by the effects of higher dimensionality 
or by the quasiparticle interaction. We model these effects below by introducing the cutoff 
energy E c < E\. 

We have found that the temperature dependence of the observed Knight shift above T c 
|H can be well fitted by (0) with 



p(x + p) oc 



1 , x < Ei 

1 + JA/E C , Ei < x < Ei + E c 

rx" ( 9 ) 

l + W — ,Ei + E c <x, 

V x — Ei 



Ei = 25meV, A = leV. 

The cutoff energy E c was chosen to be 8 meV. The corresponding function p+(x) is also 
shown at Fig. |l] by the solid line. The temperature dependence of the resulted Knight shift 
is given at Fig. || ( the solid curve (a) ). 

Second possibility for the considered system to have a small energy scale relates to the 
fact that the Og.6 compound is close enough to the metal-insulator transition ( at the oxygen 



content 06.4 )• It could result in diminishing of the DOS value at the Fermi level ||11|| . In this 
case one can expect that the small energy scale E g appears due to the strong renormalization 
effects. However, if E g is formed in a "dynamical" way, one could expect that raising of the 
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temperature up to E g should lead to the considerable changes in the quasiparticle weight 
and, consequently, in the DOS shape. This point contradicts to our initial assumption about 
essential temperature independence of p(x). 

Until now we implied the T— independent value of chemical potential. It means the grand 
canonical ensemble description of the considered system of Cu02 planes. In other words, 
the number of quasiparticles varies with temperature, which can be realized by injecting the 
carriers from the Cu-0 chains. It appears to be reasonable in view of possible important 
role of the chains in forming such one- dimensional kind of spectrum. On the other hand, 
if one assumes the number of in-plane carriers to be fixed, then the chemical potential is 
found from the equation : 



temperature shift of p or, in other words, to the effective rescaling of E\. Hence, it can be 
shown that the considered picture of the temperature variation of the Knight shift persists 
in this case, too. 

The nuclear spin-lattice relaxation rate T-f is determined by the imaginary part of the 
dynamical spin susceptibility integrated with the formfactor over the Brillouin zone. 

The INS experiments || show the enhancement of x" near the AF wave- vector Q = (71", n). 
These AF fluctuations are apparently connected with the underlying physics, though the 
detailed theory of their microscopic origin seems to be not completed yet. It is however 
widely recognized, that the different T-dependence of the relaxation rates for 63 Cu and 17 
( 89 Y) stems from the different formfactors, filtering out the AF fluctuations at O (Y) sites. 
Thus we believe, that the copper RR is mostly determined by the anomalous AF contribution 
which is beyond the scope of our paper. On the other hand, we expect the oxygen (yttrium) 
RR to be mainly determined by the ordinary Fermi-liquid contribution. We show below, 
that this latter expectation is supported by the proposed picture of DOS, essentially varying 
near the Fermi level. 




const 



(10) 



The consideration of the latter case with the above model DOS shape ([|) leads to the 
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To demonstrate it, we consider the NMR relaxation rate T{~ in the bare-loop approxi- 
mation, i.e. ignoring the effects of carriers' interaction. Then in the case of magnetic field 
along z one can write the following expression ||: 

(TTi) -1 = ^- ( dx ~^±4 - (11) 
v ; 2HJ 4Tcosh 2 (x/2T) V ; 



~p 2 (x) = E (|^(q)| 2 + |^(q)| 2 )K| 2 K +fe | 



2 



kq 



(12) 



x5{x - E k )5{x - E q+k ) 

In the case of Fermi gas and isotropic contact hyperfine interaction ( v4 a (k) = A ) one ob- 
tains p 2 (x) = 2A 2 p 2 [Ep) and the well known Korringa relation (TTx) -1 = AnTr 1 ('y n /'y e ) 2 K 2 . 
Now, when the p(x) is decomposed onto the normal and singular parts @ , one can write 

P 2 (x) ~ Pi + CxPnPsix) + C 2 p 2 (x), 

where the relative weights C\ and Ci are defined by the g-dependence of formfactors 1^4^(0.) | 2 
in flI2|). It can be seen, that the term p 2 (x) appears in (0), when both wave- vectors k and 
k+q are near the "extended" saddle point, so that the formfactors ^'^(q) have nearly zero 
argument. This observation allows us to drop the g-dependence of formfactor in evaluating 
the relaxation rate, at least at high temperatures, when p 2 plays a dominant role. 

After this remark we calculate RR as the function of temperature, taking p 2 (x) = p 2 (x). 
The logarithmic divergence at x = Ex + p, is cut off by the energy E c , discussed above. 
The obtained T— dependence of RR is shown at Fig. ^ by the dotted line (a). One can 
see from this Figure that the calculated relaxation rate is in qualitative agreement with the 
experimental findings for 17 and 89 Y 0,0. 

At last, we note, that a strongly asymmetric shape of DOS function ([|) provides a possi- 
bility to explain the almost T— independent Knight shift in the " overdoped " YBa 2 Cu 2 6 .9 5 
compound. Indeed, let us assume that the main difference of "90 K" compound from the 
"60 K" one is another value of the chemical potential or, equivalently, another value of E\ 
in®. 
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Then, taking E\ < 0, one can see from ([!]), (|TTD that i) at low temperatures the values 
of K and (TTi)^ 1 are defined by the DOS values p(fi) oc 1 + yA/\Ei\, and ii) at the higher 
temperatures the region of relatively low DOS values comes into play. 

We calculate the temperature dependence of the Knight shift and RR with ([!]), (|9]), (|Tl"D. 
similarly to the above "underdoped" case, but now taking E\ = — 25meV. The results are 
shown at Fig. ||] by the curves (b) . One sees that in this case i) raising of the temperature 
above ~ \E\\ leads to somewhat decrease of the Knight shift and RR, and ii) the limiting 
values of the Knight shift and RR at high temperatures should coincide for both over- 
and under-doped compounds, since at T ^> \Ei\ the sign of Ei is inessential. Both these 
statements are clearly supported experimentally @,||]. 

Finally, addressing to the experiments with the underdoped YBa 2 Cu 2 6 . 6 compound, 
we explain the temperature variation of both Knight shift and oxygen (yttrium) NMR 
relaxation rate within the frames of Fermi-liquid picture. We argue that the observable 
T— dependence of these values can be understood if the T— independent density of states 
in underdoped compound strongly varies with the energy near the Fermi level. ( At the 
same time the copper relaxation, apparently determined by the AF fluctuations, remains 
unexplained in our simple model. ) We show further, that the qualitative difference between 
underdoped and overdoped materials is accounted for the change of chemical potential only, 
if the shape of DOS function is delivered by the quasi-lD character of fermionic spectrum 
in YBaCuO compounds, observed in the recent photoemission experiments 
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FIGURES 

FIG. 1. The function of DOS p + (x,5) restored with regularization discussed in text (dashed 
line). The model DOS function, given by (||), is shown by the solid line. 

FIG. 2. The calculated temperature dependences of the Knight shift and (TTi) -1 with the 
model DOS function given by (^). (a) The "underdoped" case with E\ = 25meV and (b) the 
"overdoped" case with E\ = — 25meV 
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